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r^/ Introduction 



Focusing on the connection between the Landau theory of second order phase transitions 
and the holographic approach to critical phenomena, we study diverse field theories in an anti- 
de Sitter black hole background. Through simple analytical approximations, solutions to the 
equations of motion can be obtained in closed form which give rather good approximations of 
the results obtained using more involved numerical methods. The agreement we find stems from 
rather elementary considerations on perturbation of Schrodinger equations. 



^ , Much activity has been centered in the last few years on the application of the gauge/gravity 

\ duality, originally emerged from string theory [T]-[3], to analyze strongly interacting field theories 

by mapping them to classical gravity. More recently such duality was successfully applied to 
^ . describe, holographically, systems undergoing phase transitions like superconductors, superfluids 

\ and other strongly interacting systems [l]-|17j (for a wider list of references see for example |18j). 

CN ' In the holographic approach to the study of phase transitions one starts, on the gravity side, 

with a field theory in asymptotically Anti-de Sitter space-time with temperature arising from a 
black hole metric, either introduced as a background or resulting from back reaction of matter on 
the geometry. Then, using the AdS/CFT correspondence, one can study the behavior of the dual 
field theory defined on the boundary, identify order parameters and analyze the phase structure of 
the dual system. 

I Several models with scalar and gauge fields in a bulk which corresponds asymptotically to an 

AdS black hole metric have been studied [l]-|17j. The rather involved systems of non-linear coupled 
differential equations that have to be solved require in general application of numerical methods. 
In this way one often finds nontrivial hairy solutions that cease to exist for T > Tc, where T is the 
Hawking temperature associated to the black hole and Tc a critical temperature depending on the 
space-time dimensions and the parameters of the model on the bulk. 

The main point in this AdS/CFT based calculation is that the asymptotics of the solution in 
the bulk encodes the behavior of the QFT at finite temperature defined on the border. One finds in 
general a typical scenario of second order phase transition. The critical exponents can be computed 
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with a rather good precision and coincide with those obtained within the mean field approximation 
in a great variety of models. 

It is the purpose of this note to get an analytical insight complementing the numerical results, 
with a focus on the connection between the holographic approach and mean field theory for the 
calculation of critical exponents. To this end, it will be important to first stress some connections, 
already signaled in [16], between the Landau phenomenological theory of second order phase- 
transitions and the gauge/gravity approach. We will expose the importance of analyticity to explain 
the similitude in the results for relevant physical quantities (critical exponents, dependence of the 
critical temperature on the charge density and magnetic field, etc.) in very diverse models. This 
will be done first by stressing the relationship between the equations of motion and a Schrodinger 
problem, so that usual perturbative techniques allow to prove the critical behavior. We then 
show that simple matching conditions lead to results that broadly agree with elaborate numerical 
calculations. 

It should be noted that our calculations assume that the backreaction of dynamical fields is 
negligible (probe approximation) which is valid when the gauge coupling constant is large. This 
approximation is useful to study the behavior near the phase transition which is precisely the 
domain we will analyze, comparing our analytical results with those obtained numerically. In fact, 
the holographic results we have previously obtained solving numerically the Einstein- Yang-Mills- 
Higgs equations of motion both considering the probe approximation [15] and the case in which 
matter backreacts on the geometry [T7] show that when parameters are such that condensation 
takes place the solutions are very similar and the nature of the phase transition is identical. 

With this in mind and using an analytic approach proposed in [T3], we study the equations 
of motion for different models, showing how one can determine in a very simple way the critical 
behavior of the systems defined on the border with a good agreement with numerical results. 

Holography and mean field results 

In the Landau approach to second order phase transitions one considers an order parameter O and 
assumes analyticity of the free energy F = F(T, O), which can be expanded in even powers of O 

F{T, 0) = F = Fo + F2[T]0^ + F^Tp^ + ... (1) 

The dependence of the order parameter on T is obtained by minimizing as a function of O. The 
next step is to expand the coefficients F2 and -F4, in powers of T — Tc. To ensure stability and a 
change of behavior at Tc, one takes F2 = a{T — Tc) and F4 = 6/4 with a, b positive constants. In 
this way, one finds that the minimum of F is at O = for T > Tc (disordered phase) while for 
T < Tc one has 

O ~ (Tc - r)^/2 (2) 

Let us give a brief description of gauge/gravity approach to phase transitions to connect it to 
the Landau theory. On the gravity side one considers a classical field theory in a Schwarzschild-AdS 
black hole background (or one in which back reaction of fields on an asymptotically AdS space leads 
to a black hole solution) . The choice of such a geometry is dictated by the fact that the warped AdS 
geometry prevents massive charged particles to be repelled to the boundary by a charged horizon 
and as a result a condensate fioating over the horizon can be formed. To find such condensate 
one should look for non trivial static solutions for the fields outside the black hole by imposing 
appropriate boundary conditions. The behavior of fields at infinity then allows to determine the 
dependence of the order parameter on temperature, identified with the Hawking temperature of 
the black hole. 
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Basically, the free energy F is identified with the minimum of the action in the gravitational 
theory with prescribed boundary condition for the fields. Regularity of the fields at the horizon and 
smoothness of the background geometry, basic assumptions within the gauge/gravity duality, imply 
analyticity of this action. This is a first point of contact with Landau theory of second order phase 
transitions leading to a mean field behavior. We shall now see that if continuity and smoothness 
of fields are imposed in the region between the horizon and infinity, the mean field behavior found 
using a numerical approach is reproduced with a good precision, as can be seen following a very 
simple analytic approach proposed in [13] which we apply below for different models. 



The Abelian Higgs model in d space dimensions 

Dynamics of the system is governed by the action 

S = I d'^+'x^M (^-If^.F^'^ - |V^ - iA^^l^ - m^l^f^ (3) 
The background metric is the standard ^dS^+i-Schwarschild black hole 

ds'^ = -fir)dt^ + -r^dr"^ +r'^dxidx\ i = 1, 2, . . . , d - 1 

fin 

m = -^(i-^) (4) 

with 1 = 1,2,..., d— 1 and we have chosen the AdS radius to be unity. Different values for the mass 
term can be considered [6]|12j: = 0, —2 and m? = rn^p = —9/4 in d = 3 and = 0, —3, —4 
in d = 4, the last value corresponding to the Breitenlhoner-Freedman bound rn^p < marking the 
boundary of stability for a scalar field in AdS [19] . The black hole temperature is given by 

In order to look for simple classical solutions for this model one can propose the following ansatz 

^ = |M.|=^(r), ^^ = 0(r)5^o (6) 
Imposing regularity of the solution at the horizon (r = r/j) one gets 

<t^{rh) = (7) 
Concerning the asymptotic behavior of the scalar potential and the scalar field il) one has 

with 

\± = ]^{d± \/d^ + 4m2) (9) 
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According to the gauge/gravity correspondence, fi corresponds to the chemical potential in the 
dual theory defined on the boundary and p to the charge density. Concerning the scalar field tp, 
both falloffs are acceptable provided the following condition holds [20] 



- — <m^ < — - + 1 (10) 
4 4 

One can then pick either ip^ = or tp- = leaving a one parameter family of solutions so that the 
condensate of the dual operator will be given by 



{on = V. (11) 



with the boundary condition written as 



V'p = (12) 



where a, p = +, — and e"*" =1. For definiteness we shall impose the condition ip- = 
It will be convenient for the analysis that follows to change variables according to 

z = ^ (13) 
r 

so that the horizon is fixed at z = 1 and the boundary at z = 0. In terms of this new variable z, 
the equations for the functions in ansatz ([6]) become: 

- -0, (15) 



z z'^il-z^) 



where the prime denotes d/dz. 

Conditions ([7]) and ([8]) now read 



^'(1) = -(mV(i)V'(l) = 

</.(!) = (16) 
V^b(z) ~,^o D+ z^++D_ z^- 

Mz) -z-^o - qz'^-^ (17) 

with D± = V'lt and q = p/r^~^. 

In view of eq.(l5]), the system (jl4p -()15p only depends on the black- hole temperature T through 
the non-linear (fp"il) term in the scalar field equation ()14p : this describes the coupling of the scalar 
to the electric potential and gives an effective negative mass. 

In the limit T — t- oo the nonlinear term in ()14p vanishes so that the equation becomes linear with 
no non-trivial soliton solutions. Now, as already pointed out in [5] the coupling of the scalar to the 
Maxwell field is responsible for producing a negative effective mass for -0, and this effect becomes 
more important at low temperatures. This indicates that one should expect an instability taking 
place at some point towards forming scalar hair. At such point, the stabilizing effect of the m- 
term is overcome. If a non-trivial solution exist at low temperatures, with an asymptotic behavior 
associated with a non zero v.e.v. of an order parameter (O^) in the dual theory, it is natural to 
expect that a critical temperature should exist such that {0^)t>Tc = 0- This was indeed verified 
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numerically for different values of the mass m (including the Breitenlohner-Freedman bound value) 
and various space-time dimensions [6].|12j. 

Now, eq. (|14p for fixed potential <j) can be viewed as a Schrodinger equation with zero energy: 
in order to have non trivial solutions, the resulting potential has to adjust itself for the existence of 
a unique zero eigenvalue for every temperature below the critical one. In the vicinity of the critical 
point, the variation of the potential is linear in the temperature variation, but quadratic in the 
normalization of the field through the effect of the non linear term in equation (jlSp . The eigenvalue 
of the Schrodinger equation is, to first order, simply proportional to some integral of the potential 
as can be seen from its variational evaluation: with a suitable weight function a{z), the integral of 
a{z)ip{z) times the left hand side of equation ()14p gives the eigenvalue for a normalized ip and its 
variations are of second order with respect to the variations of ip. Fixing to the eigenfunction 
for T = Tc, the two sources of the variation of this integral through the potential must compensate 
themselves and one obtains a AT oc 0^ relation. It is precisely this type of behavior that has been 
found from the numerical solution to the differential equations on the gravity side. 

To go further in the analysis without resorting to a numerical analysis, we shall consider ex- 
pansions of the fields in the bulk near z = 1 and z = 0. Imposing the conditions of continuity 
and smoothness of the solutions at a point z^ intermediate between the boundary {z = 0) and the 
horizon {z = 1) will give algebraic equations between the parameters of the solution. 

For the solution near the horizon (z = 1) we have, up to order (z — 1)^: 

i^H{z) = + ^1 - 1) + ^ ''/'2 - 1)^ 

0h(^)= 00 + 01(^-1) + ^ 02(^-1)' (18) 

with ipQ, ipi, ^^2,00) 01 02 constants. The boundary conditions (fT6|) at z = 1 imply 

V'l = -(mVd)Vo (19) 
00 = (20) 

Substituting these values in (fTSj) and using the differential equations (fH|) - P^ we can obtain 02 
and -02 as a function of 0i and ipQ. We get: 



0.h(^)= 0^0 + ^^0(1-^) + ^^2^2 V'0(1-^)' (21) 



'h 

0h(^) = -01 (1 - z) + ^ (d - 3 - 2 ^^o/d) 01 (1 - zf (22) 

As announced, imposing the conditions of continuity and smoothness at an intermediate point 
Zm allows to obtain a solution. Interestingly enough, as first observed in [13] the result of this crude 
approximation are quite stable with respect to the intermediate point, so we will consider the case 
Zm = 1/2. This can be understood from the fact that ip is the ground state of a Schrodinger 
equation, so that it cannot have nodes: additional terms in the expansion of ip should rapidly fade 
away. 

We will first analyze the boundary condition Z)_ = 0, so we have the set of equations 

^^^(1/2) = ^b(1/2) , 0^^(1/2) = 0^^(1/2) (23) 
0h(1/2) = 0b(1/2) , 0'h(1/2) = 0'b(1/2) (24) 

to be solved for ipQ, 0i, D+, and /i. We obtain 
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d 16{d-2)Bq + 2'^{d-5)rhA 



2d+i 

0A+-1 



{id + 



dB' 



■"00 



B 

1 8dqB + 2'^rhA 



(for simplicity is written in terms of iJjq). Here 



A = v/l6d2 A+ + m4(A+ + 2) + 2(im2(6 + 5A+) 
S = VA+ +2 



(25) 

(26) 
(27) 
(28) 



(29) 
(30) 



Using the AdS/CFT dictionary (jlip . we can identify the v.e.v. (O) of the operator O dual to the 
scalar field with the asymptotic coefficient "0+) (O) = 4'+ = ''^^^ -^+- Now, one can write Dj^ (or 
i/'o since they are proportional) as a function of T = (d/47r) r/j. Remembering that q = p/r'^'^ one 
has 

T' 



d-l 



1 



where we have defined 



T 

24-'^(d-2) S 
{5-d){4Tr/dY-^A 



We then have, for the order parameter 



(O) = Cd 



with 



{4d + vn? 



-4m2+l)/2 



d 

2^ 



T 



{5-d) 



1/2 



(31) 
(32) 

(33) 
(34) 



(A+ + 2) 

One can see from eq. (j33p that for T close to Tc one has the typical second-order phase transition 
behavior (O) oc y^l — T /T^. Note that oc p^/^'^"^) in agreement with the change of dimensions 
of the charge density for different d's. 

Our results for the critical temperature, as inferred from (|32| ) - (|34|) in the cases d = 3, 4, m? = —2 
(with Zm = 1/2) are Tc = 0.15p^/^ and Tc = 0.2p^/^ respectively. They can be compared with those 
obtained using a different analytical approximation based on perturbation theory near the critical 
temperature. For the case m = niBF, we obtain Tc = 0.12/>^/2 a^fi j"^ = 0.25jO^'^^|21]-[22j in 
very good agreement with the exact numerical results, Tc = O.lbp^^^ and Tc = 0.25p^/^. One can 
conclude that there is a good quantitative agreement between the three sets, which is not much 
affected by the choice of the point in the method. This last fact was already observed in [13] 
for the particular case d = 3 compared with the numerical results given in [6] . 
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Non Abelian gauge field in d = 3 space dimensions 

We now consider the case of an SU{2) Yang-Mills theory in 3 + 1 dimensional Anti-de Sitter- 
Schwarzschild background, as a prototype for the gauge/gravity duality in the case of pure gauge 
theories [8]-[9]. The action is 

S = -\j d^xy^lF^.F'^'^'^ (35) 

with the field strength defined as 

F;, = d.Al - d,Al + ge'^'^A'^Al , a = 1, 2, 3 (36) 
Writing the gauge field as an isospin vector 

A = {A\ A"^ , A^) = {A^^dx^") (37) 
we consider the following ansatz for solving the equations of motion 

A = J{r)dte3 - K{r)pdipep + K{r)dpe^ (38) 
where r is the radial variable in spherical coordinates and 

X = p cos Lp 

y = p sin 99 (39) 



The background metric is 



with 



ds^ = -fir)dh + -r^dV + r'^idx^ + dy^) (40) 
f{r) 

fir) = r'(l- (41) 



In terms of the z variable defined in ()13p the equations of motion read 

t2 ^2 



r„.l«£) (42) 

1 — z-^ 

At the horizon (z = 1), J must vanish, so we have the following expansions for the fields, up to 
order {z — 1)^, 

Kh{z) = Ko + Kiiz-l) + ^K2{z- if 

JH{z) = Ji{z-l) + ^J2{z-lf (43) 

Using the equations of motion (|42p . the coefficients Ki, K2 and J2 can be written in terms of Kq 
and Ji, 



3r^ 36r^ 



JH{z) = J,{z-l)-^{z-lf (45) 
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At the z = boundary one has the asymptotic expansions 

KB{z) = Ci/rhZ (46) 
JB{z) = Do-Di/rhZ (47) 

where Di can be associated with the charge density. Coefficient Ci in ()46p should be identified 
with the order parameter {Ok) for the theory on the border. Since the order parameter is related 
to a vector field (Ai), the associated theory on the border is a p-wave superconductor [9]. 

As in the scalar field case, we will match both solutions at an intermediate point which we again 
choose as Zm = 1/2, 

KHil/2) = Kb{1/2) , K'h{1/2) = K'b{1/2) 
Jh(1/2) = Jb(1/2) , J^(l/2) = J'b{1/2) (48) 

and solve these equations for Kq, Ji, Ci, and Dq. Prom the first two identities ([I8|) . we obtain: 

4 



Jl = -^iSrl + 22Kl + 3K^lrl (50) 

(we chose Ji < so J{z) > 0). 

Substituting these values in eq. pS]) . we can obtain Kq as the root of the following polynomial 
in Kl. 

A Kl + 40 + 207rl + ASGrf, - 9rl{Dl - 48r^) = (51) 

This equation implies that vanishes as 

Dl - 48r^ 

Finally, introducing the temperature T = 3/(47r) r/j, we can write 

,2 128^2 



^0 = + 0{D\ - 48r^)2 (52) 







T 



for T near (53) 



where 

64 vr^ 

Then, for T close to Tc we have the typical second-order phase transition behavior for Ci = {Ok) oc 
^Jl — T /Tci in good agreement with refs.[8]-[9]. Our numerical value for Tc from equation ()54p . 
Tc/ V Dl = 0.091 can be compared with the numerical value obtained in [U], Td \fD\ = 0.125. 

The d = 3 scalar case in the presence of an apphed magnetic field 

We shall consider here a system with dynamics governed by the action ([3]) in the d = 3, rr? = 
—2 case, when an external magnetic field H is applied. This corresponds to the case where the 
background metric is an AdS^+i magnetically charged black hole [10j-[ll]. 

ds^ = -f(r)dt^ + -j^dr'^ + r^dxidx' , i = l,2 
f{r) 

M H"^ 

fir) = 7I-- + — (55) 
L r r 
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The /(r) = condition for having event horizons gives a quartic algebraic equation with 4 roots 
which can be exphcitly written in terms of surds. There are two complex conjugate roots which we 
shall call ri and r2 and two real roots r^, < We shall then call r/j = r4 the external black hole 
horizon. Taking L = 1 from here on, one has 



26V3 \ ' y B 

where 



r.-JrrjB.J'-^-BA ,56) 



8 31/3/^2 



The actual form of the other roots (a = 1,2,3) is not necessary since one will only need the 
standard relationships between roots and coefficients 



ri + r2 + = -rt 
rir2 + rirs + r2r3 = rl 

rir2r3 = H'^/rh (57) 



together with the relation 

Tr-73+l = (58) 



The Hawking temperature associated to the black hole is 



In terms of the variable z = Vh/r the equations of motion read 

^^ ^^ f{r,/z) "^^'^^ z^[f{r,/zy^ finjz)) 

where il^{z) = il:[rh/z) and (^{z) = (j){rh/z) (but from now on, the tilde will be omitted). In terms 
of the new variables we have: 

f{rh/z) = {rh/z)\l - z){l - nz)(l - f2z){l - hz) (61) 
df{rh/z) /dz _ -2 - Mz^/rl + 2H^z^/rl 



firh/z) z{l - z){l - fiz){l - f2z){l - fsz) 

where fa = Ta/r/i , a = 1,2, 3. 

The boundary conditions for system ()60p at the horizon are 



(62) 



= (63) 

V''(i) = ^r^wT^(i) (64) 

3rh - H^/rl 
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while asymptotically one has 

ij{z) = Diz + D2z'^ (65) 
4>(.z) = fi- {p/rh)z (66) 

For the solution of system (|60p near the horizon we have the expansions, up to order [z — 1)^: 

i^H{z) = V'O + ^1 (^ - 1) + ^ (^ - 1)^ 

cPHiz) = (Po + (l>i{z-l) + ^(l>2{z-lf (67) 



with ipQ, ipi, ip2, (po^ 01 and <j)2 constants. Using (|M|) we get 

4>o = o 



2rh 

3rh - m/rl 



V'l = ?Wt13 V'o (69) 



Substituting these values in (|67|) and using the differential equations we can obtain (j)2 and tp2 as 
functions of and ipQ. We get: 

i^{z) = + tpMtt V'o (1 - ^) - ^^(^^^^+r^'^^-^;^'^') Vo (1 - (70) 



;'i(i-^) + 



#^(1-)^ (Tl) 



As before, we impose matching conditions at z = 1/2: 

i^H{l/2) = ^b(1/2) , ^i^(l/2) = ^^(1/2) (72) 
^h{1/2) = 0b(1/2) , .^'^(1/2) = 0'b(1/2) (73) 

and solve for D2, and /i. We look for solutions with Di = and D2 corresponding to the 

order parameter of the 2 + 1 system defined on the boundary. We obtain: 

(Pi = -2R{rh,H) (74) 

_ {p/rh + 2R{rh, H)) - rj {iplm + Rivh. H) (^^g + 6)) 
^ " 2(i/2_3,4) (75) 



{3rl - H^) {p/rh-2R {rh,H)) 
2rlRirh,H) 



= A{m-3r,^f ^''^ 



where 



The equation for -f/^Q in terms of the dimensionless variable u = rh/yH takes the form 



(3n^ - 1) p/H - 2v/2/u4-6 + 7n4 

V'o = ^ , ^ (79) 

° 2u4^2/u4- 6 + 7^4 
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with the temperature ()59p given by 



y/H 47r 



1 f 1 

— 3n - ^ 



(80) 



The minimum value that u can take is the one for which T = 0, u = 3^^/^, and corresponds 
to the condition = 2^H^. From this, we see that in order to have a non-trivial solution the 

following inequality should hold 

H < — ^ (81) 
2^2/n4-6 + 7u4 

The majXiniiini of the r.h.s. is attained, for — 

> no so that there is a critical value He 
of the magnetic field beyond which no non-trivial solution exists, 

He = 0.41/3. (82) 

Using eqs. ()79p -(|80 p one can determine the critical temperature Tc as a function of H. We give 
in Figured] the resulting Tc = Tc{H) curve. Interestingly, we find that in the range 



H,>H> 



2^2/n4 -6 + 7-^4 



0.327/9 



(83) 



the curve Tc{H) becomes double valued so that a nontrivial solution only exists in the range 
> T > According to the gauge/gravity duality D2 should be identified with the order 



1.4 




0.5 



Figure 1: The phase diagram of Tc against the magnetic field H. The condensed phase {D2 7^ 0) 
corresponds to the lower left part below the line. The critical temperature decreases as the magnetic 
field grow up to the critical value He- 

parameter, {O^) = D2- We obtain the following expression 

_ 5u^-2^/p/H-2y/2/u^-6^ 
^~6u^-2 ^2(2/^4- 6 + 7^4)1/4 ^ > 
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Note that D2 becomes negative for uq < u < (2/5)^/"^ or, equivalently, for < T < 0.0316\/^- 
However, when — )• the numerator and denominator coefficients of D2 that multiply i/jq cancel 
out. This is consistent with the result of the no-magnetic field model. (In fact, we have checked 
that the results of both models, H = 0, d = 3 and iv? = —2, are identical. We present in Figure [2] 
curves for D2 as a function of temperature for different values of the external magnetic field. 




T 



0.6 0.8 1.0 1.2 



Figure 2: A plot of the order parameter D2 as a function of temperature T. The charge density is 
p = 1. The dashed line corresponds to a magnetic field H = 0.1, the dash-dot one to H = 0.2 and 
the solid one to H = 0.3. (the critical value is He = 0.6). 

The results are qualitatively in agreement with those described in [IO]-[TT] where the reported 
curves obtained by numerically solving the equations of motion are similar to those in Figure 2. 

Summary and discussion 

We have analyzed a number of models which have been proposed to study phase transition through 
the AdS/CFT correspondence. The common feature of all three models we discussed was that the 
space time bulk geometry was an Anti-de Sitter black hole. Although the dynamical field content 
was very different — a charged scalar coupled to an electric potential, the same model in an external 
magnetic field and a pure non-abelian gauge theory — the emerging scenarios are very similar and 
always include a second order phase transition with mean-field critical exponents. 

On general grounds, we were able to explain why the highly symmetric ansatze generally used, 
produce the critical behaviors seen in mean field theory or the Landau approach. Founded on basic 
principles as the connection between the equations of motion and the Schrodinger equation, we 
clarify the similarity between several relevant quantities along a variety of models. In particular we 
showed that resorting to simple matching conditions we obtain closed form solutions that signifi- 
cantly agree with the results obtained by numerically solving the exact set of equations of motion. 
This uncovers the important role played by analyticity to explain the universal behavior of certain 
physical constants. 

The method seems to work very well near the critical temperature, though it deviates from the 
numerical results as we approach T — )• 0. In this regime our approach should be refined. 

Alternative analytic calculations have been recently presented in [23] and [25] where the phase 
transition vicinity is studied solving the equations of motion in terms of a series expansion near the 
horizon. Although the approach in these works is close to the one proposed in [13] and applied here. 
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the possibility in the latter of varying the intermediate point Zm at which the matching is performed 
allows to obtain better solutions at fixed order N in the expansion, as already pointed out in |25j . 
As we have seen, in the matching approach the problem reduces to find the solution of an algebraic 
equations system and this can be done, to the order we worked here, in a straightforward way. 
Increasing the order will of course complicate the algebraic system but in view of its main features, 
it can be handled by a simple computational software like Mathematica, at least for the next few 
orders. For a large-order expansion the method followed in refs.[23] and [25] seems to be more 
appropriate. 

Although the matching method works very well near the critical temperature, it deviates from 
the numerical results as T approaches to 0. In this regime the method should be refined. In 
particular it is to be expected that taking into account the quantum fluctuations of the gravity 
theory one should be able to go beyond mean field approximation results. Also, one should consider 
generalized Lagrangians (like the Stiickelberg one considered in |14j ) leading to various types of 
phase transitions (first or second order with both mean and non-mean field behavior) as param- 
eters are changed. There is also the possibility that including fermions in the bulk model could 
substantially change the critical behavior of the theory in the bulk (see |23j and references therein) . 
The simplicity of the approach presented here, not requiring refined numerical calculations, should 
be an asset when trying to explore these more complex situations. 
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